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^ ■ Abstract 

. We report the parton distribution functions (PDFs) determined from the NNLO 

CN . QCD analysis of the world inclusive DIS data with account of the precise NNLO QCD 

corrections to the evolution equations kernel. The value of strong coupling constant 



=5 
< 
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dig (Mz) = 0.1141 ±0.0014(exp.), in fair agreement with one obtained using the earlier 
approximate NNLO kernel by van Neerven-Vogt. The intermediate bosons rates calcu- 
lated in the NNLO using obtained PDFs are in agreement to the latest Run II results. 



Account of the higher-order QCD corrections for most of high-energy processes is important 
^ ! due to the value of strong coupling constant a s is not small for a realistic kinematics. This is 
also true for the deep-inelastic lepton-nucleon scattering (DIS) process, which provides valu- 
! able information about structure of nucleon. However, since the next-to-next-to-leading order 
(NNLO) corrections were completely calculated only lately, mostly often the analysis of the 
DIS data was performed in the next-to-leading (NLO) approximation or, as the best, with the 
approximate NNLO evolution kernels derived in Ref. [1] on the basis of calculations [2, 3]. With 
the recently calculated exact expressions for the NNLO evolution kernels [4] one can improve 
available extractions of the NNLO PDFs based on the approximate evolution kernels getting 
rid of the error due to kernel uncertainty. Even so consistent extraction of the NNLO PDFs 
from the global fits including the jet production data [5] is still unfeasible since the NNLO 
coefficient functions for the jet production process are not completely calculated. In this letter 
we describe the NNLO PDFs obtained from the updated analysis of the world data on inclusive 
DIS process [6], where the NNLO coefficient functions are known and full account of the NNLO 
corrections is possible. 

We use for the analysis the charged-leptons proton/deuteron data on the DIS cross sections 
collected in the SLAC-CERN-HERA experiments [7] with the cuts Q 2 > 2.5 GeV 2 , W > 
1.8 GeV, and x < 0.75 imposed in order to reject the kinematical regions problematic for the 
perturbative QCD (pQCD) and where the nuclear corrections are particularly big. The HERA 
data with Q 2 > 300 GeV 2 were also cut off since they have marginal impact on the precision 
of PDFs obtained, but complicate the analysis due to account of the Z-boson contribution is 
required for this kinematics. The pQCD evolution input for the u-, d-, s-quarks and gluons at 
Ql = 9 GeV 2 is 

XUy(x, Qo) = l|v^ U (l - Z) 6U (1 + 72*), (1) 
* u 
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xu s (x,Q ) = -^ Vu x a °(l - x) b ™, (2) 

xd Y (x,Q ) = -^ v x a «(l-x) b «, (3) 

d 

xd s (x,Q ) = ^af(l-x) b «, (4) 

x8b(x,Q ) = ^r) s x a °(l-x)^ +b °*» 2 , (5) 

xG(x,Q )=A G x aG (l~x) bG (l+ 1 f v ^+ 1 ^x), (6) 



where indices V and 5 correspond to the valence and sea distributions correspondingly. The 
normalization factors N^ d and Aq are calculated from other parameters using the fermion 
number and momentum conservation. The value of parameter N s is defined from the condition 
that A s gives total momentum carried by the sea quarks. The value of r] s is fixed at 0.42. 
For the PDFs parameters obtained in our fit this choice provides the value of strange sea 
suppression factor equal to 0.41 at Ql = 20 GeV 2 , in agreement to the CCFR/NuTeV analysis of 
Ref. [8]. The b- and c-quarks contributions are accounted in the massive scheme with the 0(a 2 ) 
correction of Ref. [9] included. For the lowest Q/W data used in the fit the power corrections 
are important and therefore we take into account the target-mass correction by Georgi-Politzer 
[10] and the dynamical twist-4 terms in the structure functions F 2 ^ parameterized in a model- 
independent way as the piece-linear functions of x. 

Parameters of the PDFs obtained in the NNLO fit with their errors due to statistical and 
systematical uncertainties in the data are given in Table 1. These PDFs are comparable to 
ones of Ref. [6] extracted using the approximate NNLO kernel within the errors due to the 
NNLO kernel uncertainty estimated in Ref. [1]. However, at small x and Q, where the NNLO 
corrections are enhanced, impact of the new calculations is non-negligible. With the exact 
NNLO corrections the QCD evolution of the gluon distribution at small x gets weaker and 
as a result at small x/Q the gluon distribution obtained using the precise NNLO kernel is 
quite different from the approximate one. In particular, while the approximate NNLO gluon 
distribution is negative at Q 2 < 1.3 GeV 2 , the precise one remains positive even below Q 2 = 
1 GeV 2 (see Fig.l). For the NLO case the positivity of gluons at small x/Q is even worse than 
for the approximate NNLO case due to the approximate NNLO corrections dampen the gluon 
evolution at small x too, therefore the account of the NNLO corrections is crucial in this respect, 
(cf. discussion of Ref. [11]). Positivity of the PDFs is not mandatory beyond the leading order, 
however it allows probabilistic interpretation of the parton model and facilitates modeling of 
the soft processes, such as underlying events in the hadron-hadron collisions at high energies. 
The change of gluon distribution at small x/Q as compared to the fit with approximate NNLO 
evolution is rather due the change in evolution kernel than due to shift in the fitted parameters 
of PDFs. This is clear from comparison of the exact NNLO gluon distribution to one obtained 
from the approximate NNLO fit and evolved to low Q using the exact NNLO kernel (see Fig.l). 
In the vicinity of crossover of the gluon distribution to the negative values its relative change 
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Table 1: The PDFs parameters and x 2 /NDP obtained in the fit. The errors in parameters are 
obtained by propagation of the statistical and systematical errors in data. 



Valence quarks: 




a u 


0.724 ±0.027 






4.0194 ±0.050 






1.04 ±0.35 




a d 


0.775 ± 0.073 




b d 


5.15 ±0.15 


Gluon: 








a G 


-0.118 ±0.021 




b G 


9.6 ± 1.2 




7 G 


-3.83 ± 0.51 




7 2 ° 


8.4 ± 1.7 


Sea quarks: 










0.1586 ±0.0089 




a s 


-0.2094 ± 0.0044 




b s d 


5.6 ± 1.2 




Vu 


1.12 ±0.11 




b S u 


10.39 ±0.88 


X 2 /NDP 




2534/2274 



due to variation of the evolution kernel is quite big and therefore further fixation of the kernel 
at small x discussed in Ref. [12] can be substantial for the low-Q limit of PDFs. 

For the higher-mass kinematics at LHC numerical impact of the NNLO kernel update is 
not dramatic. The change in the Higgs and W/Z bosons production cross sections due to 
more precise definition of the NNLO PDFs is comparable to the errors coming from the PDFs 
uncertainties. The NNLO predictions for the longitudinal deep-inelastic-scattering structure 
function F L at x ~ 10~ 5 measured by the HI collaboration [13] also does not change too much 
since they are given by the Mellin convolution of PDFs with the coefficient functions and are 
defined by the gluon distribution at relatively big values of x. The obtained value of the strong 
coupling constant 

a s NNLO (M z ) = 0.1141 ± 0.0014 (stat ± syst), 

is in fair agreement to af NLO (M z ) = 0.1143 ±0.0014 (stat ± syst) obtained in the fit of Ref. [6] 
with the approximate NNLO kernel and to the results of the exact NNLO analysis of the 
non-singlet DIS data [14]. 

The errors in PDFs parameters given in Table 1 are calculated as a propagation of the 
experimental errors for the data points used in the fit. We calculate these errors using the 
covariance matrix method [15] taking into account statistical and systematic errors in data 
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Table 2: The NNLO inclusive rates (in nb) for the intermediate bosons production in the 
hadron-hadron collisions. The errors are due to the total PDFs uncertainties. 



pp (1.96 TeV) 
pp (14 TeV) 



W ± 
26.11 ±0.44 
197.0 ±5.3 



Z 

7.78 ±0.11 
57.7 ± 1.5 



and correlations of the latter as well. We also take into account the theoretical errors due 
to possible variations of the strange suppression factor i] s and the c-quark mass m c . For this 
purpose we re-calculate the error matrix for the PDFs parameters with ?7 S and m c released. 
Since the parameters i] s and m c are not constraint by the charged-leptons inclusive DIS data 
we confine their variation adding to the data sample two "measurements": ?7 S = 0.42 ±0.1 
and m c = 1.5 ± 0.25 GeV with the errors in these measurements representing our current 
understanding of the uncertainties in these parameters. In this approach the theoretical errors 
are included into the total error in PDFs and their correlations with other sources of the PDFs 
uncertainties are automatically taken into account. The NNLO PDFs grid for the range of 
Q 2 = 0.8 -h 2 • 10 8 GeV 2 and x = 10~ 7 1 with the total uncertainties in PDFs supplied is 
available directly 1 and through the LHAPDF library 2 . The LO and NLO PDFs grids are also 
supplied to provide a tool for checking sensitivity of different calculation to the QCD order of 
PDFs. 

The NNLO inclusive rates for the intermediate boson production at the FNAL pp collider 
and the LHC calculated using this grid and the code of Ref. [16] with corrections of Ref. [17] 
are given in Table 2. The masses and widths of the W/Z bosons were set as Mw = 80.425 GeV, 
M z = 91.188 GeV, T w = 2.124 GeV, T z = 2.495 GeV, squared sine of the Weinberg angle 
x\y(Mz) = 0.2312, squared cosine of the Cabibbo angle cq = 0.9498 [18]. The errors quoted in 
Table 2 are due to the total uncertainty in PDFs including the theoretical errors considered. 
The calculations are in agreement to the latest Run II results of Ref. [19] within the errors (see 
Fig. 2). The errors in the data of Run II are bigger than one in the calculations therefore the 
latter can be used for better calibration of the luminosity, which gives main contribution to the 
measurements error. 

In summary, we provide update of the analysis of the world DIS inclusive data on the 
proton/deuteron targets with full account of the NNLO QCD corrections including the recent 
calculations of the exact NNLO evolution kernel. The value of a s is in fair agreement to the 
earlier version of the fit based on the approximate NNLO kernels. With the exact NNLO correc- 
tions applied we observe improvement in the positivity of the gluon distributions extrapolated 
to small x and Q: Now we have gluons positive up to Q = 1 GeV, i.e. throughout kinematical 
region where the parton model is applicable. The NNLO W-VZ-bosons rates calculated using 

1 http: / / sirius.ihep.su /~alekhin / pdfa02 / 
2 http: / /durpdg. dur.ac.uk/lhapdf/ 



4 



the PDFs obtained are in agreement with the recent Run II results and can be used for better 
calibration of the Fermilab experiments in view of the uncertainty in the calculations due to 
PDFs are smaller than the experimental ones. Since these PDFs are extracted from the data 
for one single process they can be used for the quantitative studies of the PDFs universality 
that is advantage as compared to ones determined from the global fits. 
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x=0.0001 




Q 2 (GeV 2 ) 

Figure 1: The gluon distributions obtained in the different variants of the analysis (solid: the 
fit with the exact NNLO evolution; dashes: the fit with approximate NNLO evolution; dots: 
the approximate NNLO gluons evolved with the exact NNLO kernel; dashed-dots: the NLO 
fit). 
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Figure 2: The NNLO calculation of the W/Z rates for Run II at Fermilab compared to the 
data. The dotted lines give the uncertainty in calulations due to errors in PDFs; the error bars 
of the data points give the total error including one due to the luminosity uncertainty. The 
branching ratios of the W/Z leptonic decays BR(W -> Iv) = 0.107 and BR(Z -> //) = 0.034 
were applied. 
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